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The linear response to temperature changes is derived for systems with overdamped stochastic 
dynamics. Holding both in transient and steady state conditions, the results allow to compute 
nonequilibrium thermal susceptibilities from unperturbed correlation functions. These correlations 
contain a novel form of entropy flow due to temperature unbalances, next to the standard entropy 
flow of stochastic energetics and to complementary time-symmetric dynamical aspects. Our deriva¬ 
tion hinges on a time rescaling, which is a key procedure for comparing apparently incommensurable 
path weights. An interesting notion of thermal time emerges from this approach. 
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In thermodynamic equilibrium, the linear response co¬ 
efficients for perturbing forces (e.g conductivity as a 
response to an electric potential) and perturbed tem¬ 
peratures (e.g. thermal expansion coefficients or specific 
heats) may be computed with the fluctuation-dissipation 
theorem. So-called Kubo formulas relate the response to 
the unperturbed correlation between the observable and 
the entropy produced by the perturbation [T]. Out of 
equilibrium such a clear picture is lacking yet. 

For nonequilibrium systems the standard linear re¬ 
sponse to deterministic forcing has been developed via 
many approaches (see e.g. [2UT1^). In comparison, there 
is a small number of results obtained for the response to 
temperature changes [MHTT] . For example there is no 
formula to express, as a function of steady state unper¬ 
turbed correlations, the thermal expansion coefficient for 
a solid kept in a temperature gradient (e.g. as in exper¬ 
imental setups [m [19] or in models of coupled oscilla¬ 
tors [SniEI])- The construction of a general steady state 
thermodynamics needs at its heart the under¬ 

standing of temperature responses, for example in defin¬ 
ing nonequilibrium specific heats |15j . A nonequilibrium 
thermal response should as well be the subject of study 
in related fields, such as climatology [31)1132] . 

This paper introduces a theory for the linear response 
to a change of the reservoirs’ temperature. We con¬ 
sider nonequilibrium overdamped systems. Mesoscopic 
systems of this kind, including driven colloids [321 131] 
and active matter [3S], are used as paradigms in the at¬ 
tempt to generalize equilibrium concepts, such as free 
energies, within the framework of a steady state thermo¬ 
dynamics p3H33] . Moreover, they offer the possibility 
to observe experimentally the statistical fluctuations of 
energy fluxes [33113Z] • 

Our approach is inspired by a scheme based on path 
weighs and developed for the response to forces nana. 
For that theory the response turns out to be the sum of 
two terms, as in previous studies ounun]). The first one 
is half of the unperturbed correlation {OS) between ob¬ 


servable O and entropy S produced by the perturbation, 
i.e. half of a Kubo formula. The second is the correlation 
— {OK 12), where the term —K/2, of which we still have 
a less intuitive grasp, is the time-symmetric part of the 
action weighting the system’s trajectories. In order to 
avoid singularities emerging in a related temperature re¬ 
sponse [ 16 j based on a time-discretization, we introduce 
a time rescaling that significantly changes the derivation. 
As a result, the susceptibility is written as sum of unper¬ 
turbed correlations containing well-defined (stochastic) 
integrals. Moreover, an intriguing and unexpected pic¬ 
ture emerges: in S, the heat fluxes as described in the 
context of stochastic energetics [33] |351 IM] appear ac¬ 
companied by a second form of entropy production (not 
present in m), which is relevant when the system is 
coupled to reservoirs at different temperatures. 

The overdamped diffusive system we consider is de¬ 
scribed by i = degrees of freedom x = {xi}, 

evolving in the unperturbed state as 

Xi{t) = Fi{x{t))-\- ^/fi^i{t), (1) 


where every Gaussian white noise is uncorrelated with 
the others, {^i{t)^j{t')) = 2SijS{t — t'). A constant tem¬ 
perature Ti is associated to each reservoir. For simplic¬ 
ity, in particular, we choose a subset T so that Ti = T 
if i S T, which is then considered as a single heat bath 
to be perturbed. An indicator function, = 1 only if 
i G T and = 0 otherwise, is useful to keep track of the 
perturbed degrees of freedom m- We seek the linear 
response of a generic state observable 0{x) to a change 
in T, namely. 


Rorit, t') 


1 5{0{t))^ 
T 5h{t') 


( 2 ) 


where h{t) ^ 1 is the modulation of the reservoir tem¬ 
peratures, 0i(t) = Ti[\ -\- eih{t)\. Note that tempera¬ 
tures do not depend on the coordinates, hence there is 
no ambiguity in the interpretation of the stochastic equa¬ 
tion. Throughout this paper we will always consider the 
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Stratonovich convention, implying standard rules of func¬ 
tional calculus. The system may be brought far from 
equilibrium (a) by generic nonconservative forces Ti, (b) 
by different T^’s, and (c) by a relaxation from an initial 
transient condition. Indeed, given that the perturbation 
is turned on at times t > 0, the initial density of states 
Po{x) at < = 0 may coincide or not with a steady state 
density. For economy of notation we do not recall this ex¬ 
plicitly in the statistical averages, denoted by (.. .)^ and 
(...) for the perturbed [h ^ 0) and unperturbed {h = 0) 
case, respectively. 

The backward generator associated to the Markovian 
dynamics Q is 

N 

^ = with Lj = Fj {x)da,^ + Tjdl ^, (3) 

i=i 

written in a notation that will be useful later. It will turn 
out that also the following modified operator is useful to 
describe the temperature response, when the tempera¬ 
ture T of a reservoir is altered: 

L(^) ^ ^ -L, = ^ ( -F, {x)d,^ + Tdl ) , (4) 

which acts on the observables as if all temperatures were 
equal to T and all forces Fj were rescaled by T jTj. While 
L gives the derivative of a state observable O with respect 
to the kinematic time t, i.e. (LO) = dt{0), gives 
the variation of each degree of freedom with respect to 
its own thermal time Tj = tTj/T, namely, 

(5) 

such that a generalized time derivative results defined. 

Before spelling out the derivation of our results, some 
physical insights on the meaning of a temperature change 
can be gained by performing the time transformation 
Tdr = Q{t)dt in Q. Taking = 1 for simplicity, upon 
perturbation Q becomes to first order in h 

i:(T) ~ [l-/i(r)]J^(a:(r))-H VT^(r), (6) 

which clearly shows that, in such time coordinate, a tem¬ 
perature perturbation is equivalent to a force perturba¬ 
tion. The response to a small decrease in F is given by 
a theory of nonequilibrium linear response [13] as 

= -^(o{t)[x{t')F{F) - k{T')]y ( 7 ) 

where k{x{F)) = F^{x{F)) + Td^F^x^r')). If the sys¬ 
tem were in equilibrium, so as that F = —d^T-L with % the 
system’s Hamiltonian, one would expect the correlation 
functions to be invariant under a time reparametrization. 


Therefore, from Q the response to a temperature change 
in equilibrium is obtained as 

RoT{t -1') = - LH(t')]), 

which is recognized as a Kubo formula, since in equilib¬ 
rium {0{t)lFH.{t')) = [H]. Out of equi¬ 

librium instead R{t, t') depends implicitly on h through 
the time variables and no further simplification of Q 
appears possible. 

Nevertheless, the diffusive character of the system sug¬ 
gests to replace the above time change with the space 
coordinate change yi{t)/kTi = Xi {t)/v^) m , so 
that Q reads 

+ ( 8 ) 

Perturbed averages are now calculated with the path 
weights for the process y, i.e the statistical weight P^[y] 
of the trajectory {y(s) : 0 < s < t}. In particular, for all 
times 0 < t' < t and any state observable O, the linear 
response is evaluated as [42] 


RoT{t,t') = j VyO(y{t)) 


6P’^[y] 


dh{t') 


(9) 


h=0 


Here, only terms of order 0{h) are needed, hence we can 
directly linearize the Langevin equation ([^ obtaining 

m-Fi[y) + ^h[y)-^-^hyi + ^/T^^i, (lo) 


where we recognize a standard perturbing force. 


N 

fi = ^ ^ ^jUj^VjRi ~ 
i=i 


( 11 ) 


and a second one, —heiyi/2, which is atypical in that it is 
modulated by h{t). Note that in expanding the force Fi 
it is implicitly required that the values of yi are bounded. 

The path weight P^[y] is obtained with a standard 
procedure from the Gaussian path weight for 4, regarding 


^ as a functional of y via (101 


" f 1 /■* 1 

P^[y] J -Fi - 

J^ds[^y^F^ + ^{h^yJi-e^h)]'^ (12) 

(the dependence of all terms on the time s is under¬ 
stood). Plugging (12) in we encounter the mod¬ 
ulation h, which wraps in a time derivative the stan¬ 
dard result valid for deterministic perturbations, namely 
sJ^ Iodsh[y,{yi - Fi)] = -dt<[yi{yi - Fi)]. Assuming 
that the interaction forces are two body potentials, so 
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that dy-Fi = dy^Fj, by massaging the formulas and us¬ 
ing yi\h=o = Xi we finally derive 


5P^[y\ 


5h{t') 


/ i ,=0 


_ ti j 2FiXi ^ Ff 




T. 

211 . 




(t')PW. (13) 


Given our choice Ci = 1 for i G T, the response function 
of Oit) to the variation of T is thus written as 


Rorit, t') — ^ 




— I - 

47^2 


2{0{t)F,{t')x,(t')) (14a) 


^^(^) H (14b) 

_ /o(i):,,(t')L(^)F,(f')) (14c) 


-dl{0{t)xi(t')) 


(14d) 


In (14a) Ji = —FiXi is the heat flux from the *-th bath 


[38] , In addition, in (14b) there appears a novel kind of 
heat flux, 


T 




(ay>-a.) (..«) = -!:( 


( 15 ) 


^ f) 1/. 

dr' dt' ' 


which vanishes when kinematic and thermal times coin¬ 
cide, i.e. when the system is isothermal previous to the 
perturbation. The virial Vi = —XiFi of the i-th degree 
of freedom (whose average value equals Ti = T even in 
a nonequilibrium steady state [44]) is seen to act as a 
potential energy for Xj. The meaning of J™* is thus 
understood as the difference between the heat absorp¬ 
tion rate in the kinematic time and that in the thermal 
time, generated by forces dx^Vi on Xj's. Thus, the to¬ 
tal entropy flux from the selected heat bath to the sys¬ 
tem, '^i^'f Ji/T , is side by side with the entropy flux 
J™*"/T due to the heat currents installed within 
the system by the presence of different coupled tempera¬ 
ture reservoirs. These two terms are time-antisymmetric, 
that is, they change sign by going through the trajectory 
backward in time. The remaining terms, namely (14c) 


and (14d) contain the correlation between the observable 


and time symmetric quantities. 

During the last decade there was an increase of interest 
in time-symmetric fluctuating quantities (see e.g. [T^HSl - 
|49]), as it is becoming clearer that they must comple¬ 
ment entropy fluxes for a deeper understanding of statis¬ 
tical mechanics. For example, the dynamical activity of 
a jump process (counting the number of jumps) is a key 
aspect for characterizing glassy dynamics 


In linear response it was found that the time-symmetric 
sector of path weights is often related to the mean ten¬ 
dency of the system to change the perturbing poten¬ 
tial, e.g. for jump processes it becomes a shift in escape 
rates min]. The adjective “frenetic” was used to label 
this property [TU |T31 ED] . 

In (I4c) we find frenetic contributions that do depend 


explicitly on forces of the system while the last term 


(14d) does not. The presence of such term is necessary 


for having a possibly non-zero response also for free dif¬ 
fusion. For example, for a free particle starting from 
a:(0) = 0 the theory yields a response of the mean square 
displacement (a;^(t)) to an increase of T which is cor¬ 
rectly = g^df,{x‘^{t)x'^{t')) = 2 (or more in 

general twice the mobility, if we had put such constant 
in front of the forces ED; the calculation considers the 
Gaussian statistics and uses Wick’s theorem). 

Upon integration of (14) we get the susceptibility 


XOTit) = j'dt'RoTM = ^ [{0{t)S) - {0{t)K)] 

(16) 

where S, the entropy change of reservoir T, contains 


the Stratonovich integrals of (14a) and (14b), while the 
“frenesy” K contains the remaining integrals of (14c) 


and (14d|. In —K in particular there appears a pair of 
boundary terms ^ which 

left derivatives are performed to keep t' < t. 

As an example, we show the susceptibility of the en¬ 
ergy (O = R) to a change of T in a model of elastic 
slab between two thermostats, simulated using a Heun 
scheme |5S| so that the points of the trajectory can be 
used in discretized Stratonovich integrals. A scalar Xi 
is defined for i labeling a site in a portion L x L x 2 
of a cubic lattice. Each Xi in the lower L x L layer is 
thermalized at T while the Xi’s in the upper sites are 
maintained at T 2 ^ T, so that the system is out of equi¬ 
librium due to a constant heat flux. The total energy is 
T-L = J2i'>J'ixi) + uix^ - Xj) with u{x) = ^ + x 
(i -G j indicates the nearest neighbor pairs, with periodic 
boundary conditions within each layer). We compute 
XHT(t) both by direct application of a constant pertur¬ 
bation AT = T ■ 10“^ turned on at time t = 0, 

XHTV) — — - j 


and by ( |16| ). We find that the two estimates of the sus¬ 
ceptibility are in good agreement. For instance, starting 
from the system in the steady state at t = 0, Fig. [^a) 
shows the results for a slab with L = 5, when T = 2, 
T 2 = 3. Since here a positive AT brings T closer to T 2 , 
in response the mean heat current (J) from the reservoir 
T is lowered, as captured by the asymptotic anticorrela¬ 
tion between J and energy [integral of (14a) in Fig.[^a)]. 
Hence, unlike in equilibrium, the entropy flux J/T is by 
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o— oZh/ 0 ■ Eq.(16) - integral of (14a) - integral of (14c) 

—.Z^OfromAr . integral of (14b) - integral of (14d) 


o— oZ^ftO ■ Eq.(16) - integral of (14a) - integral of (14c) 

•—X^OtromAT . integral of (14b) - integral of (14d) 



FIG. 1. (Color online) Susceptibility of the internal energy to 
a change of T in the elastic slab w ith fixed T 2 ^ T, computed 
both by the direct perturbation (171 and through (161: (a) 
steady state at t = 0, and (b) transient from a configuration 
given at t = 0. Terms of (|16[) specified in (|14[| are also shown. 


for transient conditions and for steady states gener¬ 
ated by nonconservative forces or by temperature gra¬ 
dients. The understanding of the response to tempera¬ 
ture changes complements the previous approach based 
on path-weights, where the response to forces was de¬ 
rived [H [ig. We can thus say that the picture of 
the linear response for nonequilibrium overdamped sys¬ 
tems is almost complete. To fully close the circle, one 
needs the temperature response for systems with space- 
dependent noise prefactors, occurring for instance with 
hydrodynamic interactions. Investigations of these issues 
should consider a time rescaling, a key procedure in our 
approach, which leads to the concept of thermal time. 

Acknowledgments: The authors would like to thank 
K. Kroy and C. Maes for useful discussions and for their 
comments on the manuscript. 


itself not sufficient even for determining the global trend 
of the response. Fig. shows that all terms in (14) are 
relevant. To remind that the theory is applicable also 
to transient conditions, in Fig. m we show results ob¬ 
tained by starting at t = 0 from a given configuration 
with cci = 1/2 in the lower layer and Xi = —1/2 in the 
upper one. In a similar way, one might analyze data from 
a temperature quench as usually done for models of spins 
or glasses OlinillllllelllH]. 

The response formula (14) simplifies when all the reser¬ 


voirs are at the same temperature previous to the per¬ 
turbation, because = 0 and = L: 

RoT{t,t') = ^ ^{2{Omit')) + ld^Oit)x^{t')) 


ieT 


+ 


{0{t)F^{t')) - {0(t)x.(t')l.F,(t'))]. (18) 

Moreover, if the system is in a nonequilibrium steady 
state, L can be conveniently expressed in terms of the 
generator of the dynamics reversed in time, L*, as L = 
L*+2i;-Va;, with v = J / p the state velocity, i.e. the prob¬ 
ability current over the probability density [B] . Taking for 
simplicity only one degree of freedom a;, it is possible [52] 
to recast the temperature response in the form 


1 


RoT{t-t') = -—[{0{t)x{t')F{t')) 


(19) 


-b 2 {d^O{t)v(t)x{t')x(t')) F 2 { 0 {t)v{t')f{t'))]. 

In equilibrium u = 0 and F = —d^R, hence only the 

entropic terin — 7 ^ {0{t)x{f)F{t')) = 

survives, and (19) correctly reduces to a Kubo formula. 


The nonequilibrium corrections are the correlations be¬ 
tween the observable, the state velocity and the perturb¬ 
ing forces / and x. 

In conclusion, for state observables, a fluctuation- 
response relation to temperature changes has been de¬ 
rived for overdamped systems out of equilibrium, both 
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